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Q , Abstract 

^ ' It has long been known that there is a relation between boundary sine- 

. Gordon theory and thermodynamics of charge neutral two-component Coulomb 

plasma on a unit circle. On the other hand, recently it was found that open 
string worldsheet description of brane decay can be related to a sequence of 
, points of thermodynamic equilibrium of one-component plasma. Here we con- 

QJ ' sider a different decay process which is specifically described by the timelike 

' boundary sine-Gordon theory. We find time evolution to be mapped to a one- 

dimensional curve in the space of points of thermal equilibrium of a non-neutral 
^ ■ two-component Coulomb plasma. We compute the free energy of the system 

. and find that along the curve it is monotonously decreasing, defining a ther- 

, modynamic arrow of time. 



0^ 
O 
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^ '. 1 Introduction and summary 

There is an interesting proposal [1] to use the value of a tachyon field rolling down in 
^ \ its effective potential as a "clock" , when the tachyonic system is coupled to gravity, 

^ I to give a definition of time. An interesting a priori unrelated question is whether the 

concept of time could be "emergent", and associated with a large N limit. 

Previous work on tree level string worldsheet correlation functions in (open string) 
rolling tachyon background unraveled (in certain case) U{N) matrix structures [2-7]. 
It was also understood that correlation functions can be related to expectation values 
of periodic functions in an ensemble of random U{N) matrices of varying rank [5,7]. 
The ensemble can alternatively be interpreted as a grand canonical ensemble of point 
charges on a unit circle, the Dyson gas, with different points of thermal equilibrium 
labeled by different values of the chemical potential corresponding to different instants 
of time [5,8]. Later times turn out also to be related to larger values of the average 
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number N of point charges. Worldsheet correlation functions can then be related to 
thermal expectation values of moments of the electrostatic potential in the presence 
of test charges. Finally, embedding the system in a curved two-dimensional plane, 
at the critical (scale invariant point) one can reinterpret the field equations arising 
from the worldsheet beta functions. In particular, time evolution equations turn 
out to correspond to differential equations relating thermal expectation values at 
neighboring points of thermal equilibrium, at different chemical potential. 

The U{N) matrix structure was associated to a specific choice of a rolling tachyon 
background, the "half S-brane" profile T(X°) = exp(X°) on a bosonic open string 
theory worldsheet. Other choices lead to other structures. For example, generalization 
to superstrings leads to U{N) x U{M) type random matrix ensembles [6,7]. However, 
the language of random matrices may be too restrictive. For example, the random 
matrix interpretation of Dyson gas at generic (inverse) temperature was found only 
relatively recently [9]. It is easier to study what thermodynamic ensembles correspond 
to different types of tachyon condensation. 

This letter is a companion to [8] and focuses on identifying and studying the 
corresponding thermodynamic ensemble for the "full S-brane" rolling tachyon profile 
T(X°) = cosh(X°). It turns out to be the grand canonical ensemble of two-component 
plasma of equal but opposite charges confined to a unit circle on a two-dimensional 
plane. In this case there is no known random matrix interpretatioiil]. An interesting 
feature of the system is that both components of the plasma a priori have their own 
chemical potentials. Time evolution of the brane traces out a curve on the chemical 
potential plane. We evaluate the Helmholtz free energy along the curve and find it 
to be monotonously decreasing along the direction corresponding to later times. This 
"thermodynamic arrow of time" is a generalization of the one previously found in [8]; 
we also show how the latter is recovered in the appropriate scaling limit. We are left 
wondering whether there could be a more general relation: if a grand canonical statis- 
tical mechanical system with a space of points of scale invariant thermal equilibrium 
has a curve along which the free energy is monotonously decreasing, the curve turns 
out to correspond to the time evolution of some system in string theory. 

2 The full S-brane and the disk partition function 

We consider a rolling tachyon deformation of the open string worldsheet theory [10], 
called the full S-brane background. The non-trivial part of the action is given by the 
timelike boundary sine-Gordon theory (TBSG) [4] 

^0 + SSbdry = -— I dX^dX^ + \oidt cosh(X°(t)) . (1) 

27r Jdisk J 

^At least not to the authors' knowledge. 
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As was discussed in [8], the basic quantity for identifying the statistical mechanical 
system is the disk partition function (separating out the zero mode X° = x° + X'^ 
and leaving it unintegrated) , 

^disk(x°) = J M'°e-^°e-^°/^* cosh(xO(t)) _ 

Next one would expand the boundary perturbation in power series, work out the 
contractions between exponentials e^^"^*-* and sum up the series in the end. A subtle 
issue arises from the Wick ordering of the terms in the correlators, which require 
contractions of opposite types of exponentials, and cause renormalization of the cou- 
pling constant. Exploiting the underlying SU{2) current algebra structure [11] or 
proceeding via fermionization [12] shows that the bare coupling Aq gets renormalized 
teg 

A = sin(7rAo) . (3) 
After Wick ordering, all correlators are of the type 

,N+ JV_ \ 

. . . , t^,; n, . . . , r^.) ^ ( n e^'°(*'^ H e'^'^^^"^ ) (4) 

\j=l m=l / 

_ ^<i<j<N+ P ^1 lll<m<n<Af_ F I 

" a=-^in^=iie*-e-™|2 ■ 

The final form for the disk partition function reads [10] 

z,..,(.") = E E ' ' ' , ' (5) 

N+=ON-=0 iv_|..iv_. 



(2vr; 



■ ■ ■ / n n '^('^1' ■ ■ ■ ^'tN+;ri, . . . ,tn_) 

■'-^i=l m=l 



1 _ 1 



1 + Ae^" 1 + Ae-^° 
2.1 The two- component plasma 

Upon analytic continuation iX, the action ([T]) is related to the ordinary space- 

like boundary sine-Gordon theory (SBSG). More precisely, it is related to the case 
/3 = 2 of the family 

SsG,p = ^ I dXBX + \Adt cos i^Xit)] . (6) 

It has been known for a long time that the tree level partition function of ([6]) can be 
interpreted as the partition function of a two-component plasma at inverse tempera- 
ture /3, confined to a unit circle in two dimensions (see e.g. [15] for a review). 



^See also [13, 14] for additional recent discussion. 
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However, in our case we are interested in the thermodynamic system corresponding 
to the timelike theory ([T]) rather than the spacehke theory The key difference 
is that in the spacehke theory the dependence on the zero mode x = X — X' is 
oscillatory, and was trivially integrated out in previous investigations. Integration 
over the zero mode enforced charge neutrality, and [16] related the full partition 
function of the boundary sine-Gordon theory to a grand canonical partition function 
of an overall charge neutral two-component plasma, summing over the particle number 
of electrically neutral configurations only. In contrast, in the timelike theory, the zero 
mode dependence is important. Physically, in the case of D-brane decay, it reflects 
the presence of a time-dependent source; for example the disk partition function gives 
the evolution of pressure in spacetime [10]. Integration of the zero mode is associated 
with Fourier transform to the canonical conjugate variable, the total energy. 

The zero mode dependence necessitates a more general analysis of the associated 
thermodynamic system. Overall charge neutrality is no longer enforced, hence the 
grand canonical ensemble of the two-component plasma is extended to include pop- 
ulations with net electric charge. Each particle component now has its own chemical 
potential, to be related to the zero mode dependence. Since there is only one zero 
mode, it turns out to parameterize a one-dimensional curve in the chemical potential 
plane. We begin by reviewing the details of the grand canonical ensemble. 

Consider two species of particles, particles carrying positive unit charge and 
particles carrying negative unit charge, which are confined on a unit circle on 
a two-dimensional plane with positions e*** and e*"^™. They interact via the 2-body 
Coulomb potential, 

^(0i,02) = -log|e^'^^-e*^^| , (7) 

the interaction is two-dimensional while the particles are confined to one- dimensional 
motion. 

Let this system be immersed in a reservoir at inverse temperature (3 = 1/T. 
In the canonical partition function, kinetic energy contributes the usual Gaussian 
integral over the canonical momenta. We focus oijfl the non-trivial contribution to 
the partition function from the potential energy, 

N+ Af_ 



1 rw rlT + 

[ZTT) + J-w rn=l 



(27r: 



,j v^ijv_ / •■•/ Y[dtiY[dTmGf3{ti,...,tN+;Tu 

+ J-w J-^i=i rn=l 



■^The partition function factorizes, and the contribution from the kinetic energy is trivial to add 
in if necessary, hence the hterature focuses on the potential energy contribution. Further, often the 
kinetic contribution is eliminated by working in the limit of infinitely massive particles. 



4 



where 



l<i<j<N+\(^ e •'I lll<m<n<Af_ p I 



- - . (9) 

The comphcated integrals can be evaluated with the help of expansion in Jack poly- 
nomials [16]. We consider first the charge neutral case N = = A^_. The integral 
was found to be convergent for /3 < 1 so at /3 = 2 it should diverge. On the other 
hand, the grand canonical partition function 

oo 2N 

ZnentUz,P) = E Wm^^AP) (10) 

was found to have a distinct feature at /3 = 2. Previous investigations [17] had 
argued that the overall charge neutral system undergoes a phase transition from an 
insulating phase for (3 > 2 (low temperatures) to a conducting phase for (3 <2 (high 
temperatures). In the low temperature phase, the opposite charges tend to form 
dipoles and hence become insulating. As one quantitative test [17], it was argued 
that in the conducting phase, the pressure of the system is non-analytic around 2 = 
in the complex fugacity plane, while in insulating phase it is analytic around z = 0. 
Ref. [16] considered the system at high temperatures, and by a saddle point analysis, 
they found that the grand partition function has an essential singularity at (3 = 2, 



^neutral \ 



2,/5)~exp(22/(2-/3)\ _ 



However, to our knowledge a detailed understanding of the phase structure of the 
system is still lacking. 

Let us then consider the general case iV+ ^ N_. The canonical partition function 
(following [16]) is 

ZN^,N.{(^) = ZMMZNSm E iV,(iV+)iV,(iV_) ^ , (12) 

I f(A)<min(]V_|_,JV_) 

where 

Zn{(3) = , ^ ^ (13) 

rfi + f^ 



2 ' 

is the canonical partition function for a single charge Coulomb gas on a circle (the 
Dyson gas), and 

The sum in (fT2l) is over partitions with an upper bound on length and the product in 
is over the cells of a partition. Particular cases are .^7v+,o(/5) = Z^^ {(3), Zq^n^ (/3) = 
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Zn-{I3)- Consider then the special value /? = 2, where by the example (fTT!) we would 
expect the partition function (I12p to diverge. However, guided by the disk partition 
function ([S]), one can adopt a prescription to regulate the integrals, the details are 
given in Appendix. The regularized partition function simplifies miraculously to 

Zn+,n.{I3 = 2) = N+\NJ.{5N+fl + 5N.fl - 5n^+n.,o) • (15) 

Next we allow particle exchange with the reservoir and move to the grand canonical 
ensemble. Unlike in the previous investigations, we let the opposite charged particles 
have independent populations with chemical potentials so the grand partition 
function of the system is 

oo oo N- 
N+=0 N-=0 ^^ + -^^-- 



With the regulated expression (]T5l) . at inverse temperature (3 = 2 (1161) simplifies to 

Z2g{z+, z.,P = 2) = + 1 , (17) 

1 — I — z_ 

resembling the disk partition function ([5]). Note that away from the poles ai z± = 1 
the partition function is regular, in contrast to flTTl) . 

As in the previous investigation [8], relating the grand partition function of the 
thermodynamic system to the worldsheet disk partition function will require analytic 
continuation, to map chemical potential to time (the explicit zero mode of ([5])). But 
now there is a new interesting twist: the two components of the plasma come with 
two chemical potentials but there is only one physical tim^. So one must restrict to 
a subset of points of chemical equilibrium, corresponding to a one- dimensional curve 
in the two-dimensional chemical potential (or fugacity) plane. Following [8] , we start 
from f|T71) . then analytically continue to negative fugacities, which we identify with 
real values of time a;°: 

z± -z± = -Ae±^° , (18) 

which defines a one-dimensional curve z_ = X'^/z^ in the fugacity plane. 

In this way, the grand partition function f|T7|) is mapped to the disk partition 
function ([5]). Note that we need to restrict to values \z±\ < 1 which correspond to 

< ^ = -logA . (19) 

It is interesting that the real time range Ax^, corresponding to the allowed values of 
the fugacities, has a very natural physical interpretation: the lifetime of the bran^ 



■*We thank S. Brodsky for this remark. 

^/.e., the time interval during which the unstable brane first forms and then decays, see e.g. [18]. 
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2.2 Thermodynamic arrow of time 

In the case of the half S-brane, it was found that iV, the average number of charges 
in the Dyson gas, could be used as the time variable, and the free energy of the 
system was a monotonically decreasing function of N, giving an arrow of time. This 
analysis was done at the temperature (5 — 2 corresponding to the exactly marginal 
rolling tachyon deformation [8] . Here wc extend the calculations for the full S-brane. 
Extension to superstrings will be considered in [19]. 

The average particle numbers for positive and negative charges are 

Restricting to the curve Zj^Z- = A^, the numbers N± are related by 

= ■ (21) 

To specify the thermodynamics of the system onto the curve, we solve the constraint 
and express the grand partition function in terms of 

(1 _ A2)e^''+ 

Then, on the curve the Legendre transform to the Helmholtz free energy is given by 

Aurve(A^+) = -P-' log Z2G + Z^" log Z2G • (23) 

At P — 2, the free energy along the curve becomes 

Acurve(A^+) = ^logiV+ - ^ log(l + (1 - X')N+) - ^7V+log(l + (1 - A^)7V+) 

+^iV+log(A^ + (1 - A^)iV+) - ^(YT:|^log(A^ + (1 - A^)iV+) . (24) 

It is straightforward to check that the free energy Acmve is indeed monotonously 
decreasing as a function of the curve parameter 7V+. So it can be interpreted as a 
thermodynamic arrow of time, in the sense of [8] . 

2.3 Dyson gas as a scaling limit of the two-component plasma 

Recall that to get half S-brane, we are shifting the origin of the time coordinate as 
follows 

Ao J dt cosh(X°(t) + C) = ^Ao J cit(e^°+^ + e"^"-^) dte""" , (25) 



where the hmit 



Ao , C ^ oo , Aoe'^ =2X= fixed (26) 

was apphed. This hmit is also reahzed in Dyson gas. RecaU the partition function 
f fT7|) . First we implement the shift —>■ + C: 



Ae±(^"+^) . (27) 



and use a different notation for the renormalized coupling, A = sin(7rAo). Now we 
apply the limit (126|) for the two-component plasma, 

A ^ ttAo ^ (28) 
z_ = Ae-^"-^ ^ (29) 
= Xe''+^ vrAoe^e"" = 27rAe"" = z , (30) 

so the population of negative charges goes to zero while the positive charges remain. 
In more detail, 

iV+ ^ . 5N+^^ (31) 

I- z ^Jz 

while 

iV_ ^ ; (5Ar_ ^ . (32) 

Also the Helmholtz free energy fl2^ reduces to the previous result obtained in the 
case of the half S-brane 

v4eurve(iV+) ^ -^[(iV+ + l) log(iV+ + l) - iV+ log iV+] . (33) 



(34) 



In particular, the grand partition function reduces to 

^2G = Zg , 

1 — z 

the partition function of the one-component plasma (Dyson gas). 
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A Partition Function at /3 = 2 

The purpose of this appendix is to fill in the gaps between (fT6l) and (ITTj) . Let us 
begin from (fT6|) . It is straightforward to manipulate 

Af+=OAf_=0 Jv + .iv_. Ar+=lAr_=l iv + .iv_. 

+ 1 + -. 1 , (35) 

1 — z+ I — Z- 



where 



2 /•TT /•TT 

Zn+,N^ = xjy , jy_ / ■ ■ ■ / dtl - ■ ■ dt^^dTl ■ ■ ■ dT]\f_ 

{ZTT) + J —IT J —IT 



(27r 

nl ZD^^Z yD'^^i I 2 I I I /D^'J"m yD^'^n I 2 

l<j<j<Ar+ p I lll<m<n<Ar 



a=^in^=iie^*— e--|2 ■ ^^^^ 

The only non-trivial part is to show the vanishing of Z]\f^^N_ with N^,N^ ^ 0. 
To achieve this, define Zi = e**' and Wm = e*"^". Then, (l36l) becomes 

(— / / dzN+ dwi dwN_ 



(27r)^++^- Ic Jc zi zn+ wi wn_ 



(37) 



^ ■ ■ ■ ^ (i^i ■ ■ ■ dzN+dwi ■ ■ ■ dwN_ 

N+-N- I ^ \N--N+ 



\Zi---Zn+) \Wi---Wn^J 

ni<i<j<Af+ ~ ^j) ] Ill<m<n<N- [~(^m — Wn) ] 



(3^ 



where C is the unit circle (integrated counterclockwise). 

Due to the double poles, one needs to properly define a regularization prescription. 
We use the principal value prescription, defined as 

11111 1 ^ 

+ o77^ ^ = P-, > C = l + e. (39) 



{Zi - Wm)'^ 2 {Zi - CWmY 2 {C Zi - WmY (Zi - Wm)'^ 

Let us now integrate over one variable, and without any loss of generality we can 
choose zi, assuming < iV_|^. Now 

zn+,n. ~ / dzi n [-(2:1 - Zjf] ■ n 

Zi j=2 m=l . 



-P- ' 



[Zi - WmY _ 

f^dz,g{z^) = J . (41) 



^If > N- we will choose wi. 
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At infinity 

hence there is no residue at 2:1 = oo. At zero there is no residue as well, since 

giz,) z^'-''^ . (43) 

We thus conclude that all the residues of g are lying on the integration contour. The 
principal value prescription gives one half of each residue contribution, 

^Af+>o,Af->o ~ J = Tri ^ Res g , (44) 

where 

5^ Res ^ = (45) 
as the sum of the residues of any rational function is zero. 
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